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A new type of transverse instability in dispersively nonlinear optical cavities, called the optical 
whistle, is discussed. This instability occurs in the mean field, soliton forming limit when the cavity 
is driven with a finite width Gaussian beam, and gives rise to oscillation, period doubling, and 
chaos. It is also seen that bistability is strongly affected due to the oscillation within the upper 
transmission branch. The phenomenon is interpreted as a mode mismatch in the soliton formation 
process and is believed to have broad applicability. 



INTRODUCTION 

Over the last 25 years the nonlinear optical resonator has been a rich source of interesting phenomena These 
systems are attracting attention recently due to their potential applications as all-optical switches, memories and 
logic gates j|, at both the classical and quantum levels. The resonator geometry naturally gives appreciably stronger 
nonlinear effects for a given incident beam intensity than do the travelling wave schemes, potentially allowing the use 
of faster, less lossy, materials with relatively weaker nonlinearities (for example, atomic vapors or even silica as opposed 
to thermal or photorefractive materials). The resonator geometry also leads to a variety of time-dependent behaviors 
which must be adequately characterized, either to avoid undesired effects in engineered systems or to perhaps take 
advantage of them for communications or computing purposes. 

The most fundamentalphysical process in the nonlinear resonator, plane wave optical bistability, was first observed 
experimentally in 1975 p}. Ikeda H later showed in a discrete return-map analysis that as the incident intensity 
increases the intracavity field intensity can oscillate and undergo a period doubling route to chaos. In the 1980's a 
great deal of work concentrated on transverse structure and dynamics in nonlinear cavities, in order to incorporate 
the phenomena of self-(de)focusing and diffraction. McLaughlin et al M demonstrated in 1985 that a finite-width 
Gaussian beam incident on a cavity with a single transverse dimension can give rise to a transverse oscillation in the 
cavity field, essentially due to a modulational instability. These oscillations are of period equal to a multiple of the 
cavity round-trip time and thus correspond to the occupation of multiple resonator longitudinal modes. More recently 
Haelterman and Vitrant discussed a transverse oscillation which arises when a nonlinear Fabry-Perot is illuminated 
obliquely, due to the lateral drift of the field within the cavity @ . 

Here we report on the discovery of a new type of oscillation within high-finesse nonlinear cavities driven by Gaussian 
beams. This optical whistle oscillation occurs as a competition between spatial soliton formation || M and the driving 
input, and can exhibit oscillatory or chaotic behavior. 

GENERAL CONSIDERATIONS 

The dynamics of a general (dispersively) nonlinear optical cavity driven with a finite-width beam is principally 
governed by four effects: propagation, damping, diffraction, and nonlinearity. A qualitative understanding of cavity 
behavior is gained by simply comparing the characteristic timescales of these processes, given by: 

Tprop L/C 
Tdamp LJ~ fC 

Tdiff = w 2 /Xc 
T nt = X/cAn 

where L and T are the cavity length and finesse, A is the optical wavelength, w is a characteristic size of transverse 
features in the optical beam, and An = \n(I) — n(0)\ is the typical nonlinear index shift caused by the optical field of 
typical intensity / (the linear index is assumed to be unity). 
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The timescale T n i is the time required for the nonlinearity to induce a 2tt phase shift in the light field. For self- 
focusing nonlinearities we can also identify t„i as the typical formation time of transverse structure arising from the 



modulational instability 10|. These structures form with characteristic transverse wavelength A = A/v An. A feature 
of this size will diffract apart in a time A 2 /(Ac) rs T n i, and we generally conclude that these structures spontaneously 
form so that diffraction and nonlinearity arc equally important. The same result applies to the formation of spatial 
solitons, which can be considered modes of self- induced waveguides |jll[ (f2). Soliton formation is more efficient when 
the incident beam is of size A, in analogy with mode-matching in linear optics. The optical whistle instability arises 
from a "mode mismatch" between the incident beam size and the (self-consistent) soliton dimension A. 

Returning to the cavity timescales of Eq. [I], one important case is when T prop is much shorter than the others. In 
this mean field limit the field changes little in traversing the cavity once, and we can suppress the field envelope's 
longitudinal dependence, ft can be shown that in this limit only a single longitudinal cavity mode is appreciably 
occupied |f~3| ]. 

The optical whistle phenomenon exists under the same conditions as spatial soliton formation: 

Tprop Tdiff ~ Tnl ^damp (2) 

For the remainder of the discussion we confine ourselves to this case. From Eq. [j] it follows that we are discussing a 
high-finesse cavity. 



THE NONLINEAR CAVITY EQUATION 

We now restrict our attention to a Kerr nonlinear cavity having a single transverse degree of freedom x. The cavity's 
internal field envelope is governed by the Lugiato-Lefever equation , written here as 

Bf ir B 2 f 

~dt=2k^~ T ' £ + 1{AU ' )£ + T ' £d " ve + iujn2 \ £ \ 2£ ' ( 3 ) 

where E is the internal cavity field envelope amplitude, k = 2ir/X is the longitudinal wavenumber, V = cT/L is the 
amplitude decay rate from the cavity (T is the amplitude transmission coefficient at each mirror, assumed equal, and 
L is the cavity length), u> is the field angular frequency, Aw' = lu — cu cav is the detuning of the driving field from 
linear cavity resonance, and 712 is the nonlinear index inside the cavity. A time-independent version of this equation 
was also derived by Haelterman et al from their modal theory [HJ . Equation || is slightly different from the version 
in Ref. H], which has been rescaled to yield V = 1; here we will employ a different rescaling to easily accomodate 
the important limit T — > 0. 

Equation |^ is now made dimensionless by choosing an arbitrary distance scale xq and relating the time and field 
scales to and |£o| to it using to — kx^/c and \£q\ = f / /exo y^jnjf. After rescaling Eq. || takes on the dimensionless form 

* = ^V 2 <F + ir]^^ + iAuV - r(* - m drlve ) , (4) 

where r\ = +1 (— 1) for self-focusing (-defocusing), "F = £/|£o|> Aw = Aw'to, and T = T'to. The conditions of Eq. || 
have allowed us to use the mean field limit in deriving this equation and also imply that T <C f . 



SOLITON FILTERING AND THE OPTICAL WHISTLE 

Equation ^ is a damped, driven version of the nonlinear Schrodinger equation. In the limit T — > and for rj = +1 
it admits the family of stationary soliton solutions [[l6| 

V(X,T) = d^sechiX/d) , . 

Alu = -l/(2d 2 ) (5J 

where d is an arbitrary width and X — x/xo, T — t/to are dimensionless coordinates. There are other, time-dependent 
("breather"), soliton solutions [ fL7| which are unimportant for our present purposes. Equation ^| is also a solution in 
the damped case (T ^ 0) if the driving field has the special form 

*drive,sech{X, T) = T) = cT 1 S6ch(A/d) . (6) 
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Whether this T ^ solution is a soliton in the mathematical sense is a subtle issue and not resolved here. We will 
take the common practical approach and refer to this as a "soliton solution" in reference to its observed soliton-like 
robustness, keeping in mind the potential inappropriateness of this terminology. 

A general qualitative feature of solitons is their robustness to small perturbations. Hence we might expect that the 
Gaussian driving field 

■$ drl ve,Gau SSl an{X : T) = d- 1 exp[-X 2 /(1.699d) 2 ] . (7) 

may lead to the soliton solution (||) nearly as efficiently as the sech drive (^) when r< 1. The constant 1.699 has 
been chosen to minimize the difference integral 

drive ,Gaussian 

{X)-y drive>sech (x)¥ dx. (8) 

This possibility was considered in a travelling wave geometry by Burak and Nasalski |l8f| , who found that as much as 
99.5% of the incident power is converted into a soliton. 

In our present situation it is found that the folded path arising from the cavity configuration can modify this picture 
substantially. Generally, soliton forming behavior is observed when the driving field amplitude is matched to its width, 
as in Eq. 0. Because the non-sech incident beam profile yields a sech beam in transmission, we call this process soliton 
filtering in analogy with ordinary spatial filtering. A larger driving amplitude leads to a "mode mismatch" since the 
soliton is trying to form with a smaller width than the input beam. The optical whistle oscillation results from a 
competition between the soliton trying to shrink in size, and the driving field feeding in light mismatched in phase 
and spatial profile. 

Considering again the Gaussian driving field of Eq. [?], recall that the choice of distance scale Xq in deriving Eq. |^ is 
arbitrary, so we can choose it to be the waist size of the input beam without loss of generality. (In this case the time 
scale to is the time of flight through the beam's Rayleigh range.) Then we have ^ d rwe(X, T) — q exp (—X 2 ) where q 
is a real, dimcnsionless driving amplitude, and our model is characterized by the four real dimensionless quantities 77, 
r, Aw, and q. We expect optimal coupling into the soliton when 



q w 1.699 
,w w -(1.6! 

The former is remarkably close to the optimal value 1.69 in the travelling wave case, found using the inverse scattering 



Aw w -(1.699) 2 /2 = -1.443. ^ 



transform |18|. 



NUMERICAL RESULTS 

Equation [| is solved using a version of the popular split-step technique |l^] [^0| , accurate through two orders in the 
time step size h. The results presented below are calculated using h — 0.01 and a spatial grid consisting of 256 points 
separated by 0.156 dimensionless units. These values are found to give accurate results for the cases of interest here; 
halving the grid spacings in space and time yields essentially identical results for a variety of test cases. 



Oscillation 



Figure |l] summarizes the asymptotic (long-time) cavity solutions when the Gaussian driving field is suddenly turned 
on at t — 0. We have here assumed an initially empty cavity, n = +1, and T = 0.14. The asymptotic cavity solutions 
are categorized into five basic types: (1) steady state solutions, (2) period 1 oscillations, (3) period 2 oscillations, 
(4) long period oscillations, and (5) chaos. The regions are seen to have complex boundaries, and points near the 
borders of the "chaotic" regions show particularly interesting behaviors: at the borders with steady state regions lie 
oscillations of very long period, and at the borders with the normal oscillating regions period doubling occurs. Period 
3 oscillations are also observed. 
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FIG. 1. Asymptotic behavior of cavity field, for rj = +1, T = 0.14, and an initially empty cavity. Labels indicate regions of 
steady state, period 1 oscillation (1), period 2 oscillation (2), and chaos (C). 

Time series of examples are shown in Figure |2|, where we plot the total power output for a variety of cases. For 
the period 1 example, snapshots of the cavity field amplitude are shown in Figure |. The field profile narrows as the 
soliton forms, then is attenuated as the nonlinearity shifts the peak out of cavity resonance and power enters at the 
edges. The fundamental oscillation frequencies for several cases are plotted in Figure |^. It is not known what ranges 
of r permit oscillations to occur. 
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(a) Acq = -0.6, q = 4.0 



(b) Acq = -0.2, q = 4.0 
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(c) Acq = -0.32, q = 4.2 
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(d) Acq = -0.32, q = 6.0 
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FIG. 2. Time series of transmitted power for specified parameters, illustrating the asymptotic solutions shown in Fig. 
period 1, (b) period 2, (c) long period, and (d) chaotic oscillation. r\ — +1, V — 0.14. 




FIG. 3. Cavity field amplitude profiles for the period 1 oscillation shown in Figure |2|( a). Alo — —0.6, q = 4.0,r = 0. 
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FIG. 4. Fundamental frequency of oscillation as a function of driving amplitude q. The solid curve is for T — 0.14, Aui = —0.2 
and shows one of the period-doubled regions of Fig. |l| The long dashed curve is for T = 0.14, Au — —0.6. The short dashed 
curve is for V = 0.3, Acu — —0.6. 



Bistability 



An interesting issue is the effect the optical whistle has on dispersive optical bistability. Figure g compares a 
bistability curve from our time-dependent analysis with that given in Ref. EJ which assumed a steady state field. 
Our respective equations have been made dimensionless in slightly different manners; to correspond to the Xq = 1, 
A = +3 curve in their Figure 3a we have chosen T = Xq/2 = 0.5, Aui = — TA = — 1.5 and have related their field 
amplitude / to ours using * = VTf. Part of the upper branch of the bistability curve is seen to be unstable due to 
the whistle oscillation, and the oscillation amplitude is indicated by the bars in Fig. |[ In general we observe that 
where bistability occurs (Aw < —0.19 for T = 0.14) the whistle oscillation occurs only in the upper branch of the 
bistability curve. However, the oscillation can also occur when there is no bistability, as is evident from Fig. [P. 
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FIG. 5. Comparison with steady state theory. The dashed bistability curve is from Fig. 3a of Vitrant et al (1990) and the 
solid curve is from our time-dependent analysis. Bars on the solid curve show the oscillation amplitude. T = 0.5, Auj = —1.5. 

In other cases the optical whistle modifies bistability curves more dramatically. Figure ^| shows steady state and 
time-dependent bistability curves for T — 0.14, Aui — —1.443. In this case the whistle oscillation has moved the 
switch-off point substantially, a common occurrence. In contrast, the switch-on point has not been seen to change 
from its steady state value in our analysis. 
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FIG. 6. Comparison of steady state (dashed) and time-dependent (solid, with oscillation amplitude indicated by bars) 
analyses. In this case the whistle instability changes the switch-off point in the bistability curve. F = 0.14, Aui — —1.443. 
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Soliton Formation 



We argued above that the soliton solution of Eq. g would be stable when r -C 1, in the sense that a sech-profilcd 
cavity field would result when a Gaussian driving beam is applied. This is numerically demonstrated in Figure 
where steady state field solutions are presented for the cases T — 0.1 and T = 10. The former solution is very close 
to the sech form of Eq. ||, and the latter is close to the Gaussian driving field. Although the existence of solitons in 
a damped cavity is an unresolved technical question, the sech solutions are seen to be quite stable when f < 1. 




X 

FIG. 7. Steady state field solution with Gaussian drive, plotted on linear and logarithmic scales, in the upper transmission 
branch. The dashed curves are the driving function and 1.58 sech(1.575x); overlapping solid curves are solutions with V = 10 
and r = 0.1, respectively. The latter solution is very close to the soliton form of Eq. ^ despite being driven by a Gaussian 
beam, providing evidence for the stability of the sech solution when r<l. q = 1.8, Au) = —1.2. 

In the travelling wave case soliton formation with Gaussian driving beams can be very efficient; as much as 99.5% 
of the input power can be transformed into a soliton |l8( ] . In a cavity we might expect this to be an upper bound since 
imperfect interference (and reflection of the incident beam) would seem to present an additional loss of efficiency. In 
Eq. |^ we made a simple estimate of the optimal coupling parameters. However, Fig. ^ shows that for Auj — —1.443, 
r = 0.14 the whistle disrupts the bistability curve's upper branch near q = 1.699. A nearby case (Aui = —1.2, 
r = 0.14) without the disruption is shown in Fig. ||. At the switch-off point q as 1.6 the coupling efficiency into the 
transmitted beam is 97.6%, still quite high. Because T — 0.14 <C 1 the transmitted field is found to be quite close 
to the sech form of Eq. ||; when T > 1 the transmitted field is more Gaussian in profile. Figure || shows the steady 
state soliton width in the upper branch of Fig. ^ and the decrease in width as the oscillation point is reached. For 
the present work we have not attempted a thorough study to optimize soliton coupling efficiency but wish to point 
out that fairly high efficiencies are possible in cavities. 
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FIG. 8. Bistability curve showing nearly complete (97.6%) coupling into the transmitted sech beam. The expected optimal 
case, Eq. 0, is disrupted by the whistle as shown in Fig. pi Alu = —1.2, T = 0.14. 
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FIG. 9. Width d of soliton solution in upper branch of Fig. ||as a function of driving amplitude q. As q increases the soliton 
narrows in width and begins oscillating at q ~ 2.4. Auj — —1.2, V = 0.14. 



DISCUSSION 

A variety of instability phenomena similar to the optical whistle have been reported in the literature [0 . These 
can be broadly categorized into multi-mode and single-mode (mean field) effects. In the former case the plane wave 
oscillation phenomenon has been called "self-pulsing" (it occurs in both the context of absorptive |2^] and dispersive 
H J2j| bistability). The instability here is essentially due to the interplay of two different timescales, nonlinearity 
and cavity feedback, and gives rise to an oscillation period which is a multiple of the cavity round-trip time. Slightly 
later the single-mode case was also shown to yield oscillations. Ikeda and Akimoto |^4| first demonstrated this within 
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a plane wave, purely dispersive Kerr model. Lugiato et al |25J] constructed a more realistic optical Bloch model and 
found oscillations (and chaos) in that case as well, again at the plane wave level of description. 

Our work is more closely related to the transverse instability discovered by McLaughlin et al p(J and discussed 
in the mean field limit by Lugiato and Lefever In these studies, stationary transverse structure in a broad 

illuminating beam arises from a modulational instability. McLaughlin et al Q showed in 1985 that the transverse 
features can oscillate in time. In all cases, more recent numerical work has shown that these transverse perturbations 
can grow to become soliton chains |?J p7| which may or may not be stationary in time. Lastly, for a plane wave 
driving field of variable intensity a series of bifurcations similar to the ones seen here are observed to occur in the 
mean field case, eventually leading to spatio-temporal chaos [ p8| . Our present work is different from these earlier 
results because here the most unstable transverse wavelength is comparable to the driving field width; the instability 
consists of a broad spectrum of transverse wavenumbers and it is somewhat more intuitive to view the process as a 
mode mismatch in the soliton formation process rather than as a sinusoidal modulation in a very broad beam growing 
into a (perhaps) time-varying soliton chain. 

There are several possible experimental configurations in which the optical whistle might be observed. One concrete 
realization is a fiber-loop geometry p9| , in which Gaussian pulses from a laser are injected into a (nonlinear) fiber 
loop and the copropagating longitudinal coordinate z — vt plays the role of the tranverse dimension x in our analysis. 
Another possibility is a cylindrical nonlinear cavity with strong single-moding in the y and z dimensions |30f . 

The optical whistle phenomenon presented here is believed to exist beyond the Kerr nonlinear, one dimensional 
model, whenever one is in the soliton forming limit of Eq. ^. We conjecture this on the basis of our general picture 
that the whistle is the result of a mode mismatch between the generated soliton and driving field. The study of related 
phenomena in other systems would be an interesting topic of further research. 
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